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ABSTRACT 



Recent improvements in rare earth magnets have made it 
possible to construct strong, lightweight, high horsepower 
DC motors. This has occasioned a reassessment of electrome- 
chanical actuators as alternatives to comparable pneumatic 
and hydraulic systems for use in flight control actuators 
for tactical missiles. This thesis develops a lcw-crder 
mathematical model fcr the simulation and analysis of brush- 
less DC motor performance. The model is implemented in CSMP 
language. It is used to predict such motor performance 
curves as speed, current and power versus torque. 
Electronic commutation based on Hall effect sensor posi- 
tional feedback is simulated. Steady state motor behavior 
is studied under both constant and variable air gap flux 
conditions. The variable flux takes two different forms. 
In the first case, the flux is varied as a simple sinusoid. 
In the second case, the flux is varied as the sum of a sinu- 
soid and one of its harmonics. 
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I. INTRODUCTION 



Direct current electric motors are gaining wider use in 
space applications where long life, high reliability, high 
torgue and light weight are critical factors. In recent 
years, advances in magnet technology and materials, espe- 
cially the use of rare earth magnets, have made possible 
significant improvements in the torgue to inertia ratios of 
permanent magnet (PM) motors. Samarium cobalt is one such 
rare earth magnetic material. The great potential of the 
rare earth magnets derive from their inherent high flux 
density and high coercivity properties. Higher flux densi- 
ties mean greater developed motor torgue while high coer- 
civity means greater innate resistance to demagnetization 
which permits thinner magnet sections. The two factors in 
combination result in increased mechanical power and energy 
product with decreased physical size and weight. 

Brushless dc motors belong to this class of permanent 
magnet motors and thus enjoy their improved torgue to size 
characteristics. In fact, brushless dc motors are signifi- 
cantly smaller than either their ac counterparts or conven- 
tional brush-type dc motors on an equal horsepower basis. 
Compared to its conventional counterpart, a brushless dc 
motor’s structure is "inside out." In the brushless config- 
uration, the permanent magnets are attached to the rotor and 
the conducting coils are placed in the stator. As its name 
suggests, brushes have been eliminated from this type of 
motor. Instead, commutation of current in the stationary 
armature is accomplished electronically by switching cn the 
appropriate windings via transistors specially designed to 
provide the high currents needed for motor power. In order 
to seguence the current to the coils, the commutation 
circuitry is provided rotor position feedback from sensing 
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devices, most commonly from Hall effect sensors. These are 
simple semiconductor devices that develop a polarized 
voltage depending upon the magnetic field passing through 
them. Small, highly sensitive and reliable. Hall effect 
sensors require little power to operate and thus are widely 
used . 

Erushless dc motors possess all the advantages tradi- 
tionally associated with conventional dc motors such as 
linear torque-speed characteristics, excellent response 
times and high efficiency. In addition, there are a number 
of important advantages that brushless dc motors possess 
over standard brush- type dc motors. An obvious advantage is 
that electronic commutation means the elimination of commu- 
tators, the wear and tear of brushes, and the build up of 
carbon dust which heretofore significantly limited the life 
of the standard dc motor. Since electronic current 
switching essentially involves no component wear, there 
results improved reliability and enhanced lifespan. As well 
as being subject to wear and tear, the action of brushes 
sliding over commutators also generates radio frequency 
interference (RFI) which is a severe liability in many 
applications. In fact, it is often in hostile and explosive 
environments that the advantages of dc motors are most 
needed. Another benefit is in heat transfer characteris- 
tics. By placing windings in the stator the thermal path to 
the environment is made shorter and more direct. With the 
removal of heat thus enhanced, mechanical and electrical 
malfunctions due to heat are reduced. A further important 
advantage of the brushless dc motor is in the area of motor 
speed variability. Whereas ac motors usually operate at one 
speed fixed by the powerline frequency, any dc motor’s speed 
can be changed to meet varying power requirements by simply 
adjusting the dc input voltage. The electronic commutation 
scheme cf the brushless dc motor lends itself nicely to just 
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this type of variable speed control. I t is a relatively 
simple matter to attach to the commutation circuit beard a 
small, inexpensive microprocessor whose logic can be 
programmed to run the motor at any of several speeds or in 
either direction. A single motor can thus be programmed, in 
place if desired, sc that its speed, direction and tergue 
are matched to a particular application. 

There are at least as many applications for brushless 
dc motors as there are for brush-type motor systems as well 
as a hest of new applications. Of particular relevance to 
this thesis is the use of brushless dc motors in flight 
control actuator sytems of Navy tactical missiles, in 
particular, the cruise missile. There are three general 
types of actuators that are in use in missile systems: pneu- 
matic, hydraulic and electromechanical. The maneuvering 
requirements of the mere advanced tactical missiles call for 
high control torques which heretofore could only be produced 
by high pressure pneumatic or hydraulic actuators. Beth of 
these actuator systems are more susceptible to mechanical 
malfunction than their electrical counterparts. Kith the 
improved torque to inertia ratios resulting from better rare 
earth magnet materials and the higher mecnanical reliability 
due to removal of brushes, as well as other advantages 
previously enumerated, the brushless dc motor is in a posi- 
tion to compete with most pneumatic and hydraulic systems. 
This thesis is but one part of a detailed study of the state 
of the art in electromechanical actuators as applied tc the 
cruise missile system [Ref. 1], 
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II. NATURE OF THE PROBLEM 
A. PROGRAMMING LANGUAGE 

The objective of this thesis is to create a simplified 
low-order mathematical model that will accurately simulate 
the- response of a brushless dc motor. The IBM Continuous 
System Modeling Program (CSMP) was chosen for use since it 
is a user-oriented computer language specifically created 
for modeling dynamic physical systems and is one of the most 
widely used in this area £Ref. 2]. It has the flexibility 
of determining the response of a system that is modeled 
either in a block diagram format or as a set of ordinary 
differential equations. While CSMP is an offshoot of IBM’s 
Digital Simulation language, it embodies Fortran as its 
source language and retains much of Fortran’s capability and 
flexibility. CSMP is an applications-oriented prcgrara 
intended for use with the IBM mainframe systems. The 
utility and ease of using this program language stems from 
simplified program control statements that almost exactly 
describe the mathematical equations or physical variables of 
the system and from the flexibility of its program structure 
which contains preprogrammed function blocks that obviate 
the need for complicated subroutine programming. In 
essence, CSMP was chosen because it permits the user to 
concentrate on the details of the physical system rather 
than on the time-consuming complexities of programming. 



E. SISTEM BLOCK DIAGRAM 

The block diagram format was chosen as the basis for 
the CSMP program since it is a convenient way of analyzing 
the input-output relationships of a physical system and the 



12 



flow of signals through it without having to refer to the 
system eguations. The basic idea of a block diagram steins 
from the application of Laplace transforms to a system's 
integro-diff erential cr differential equations. In effect, 
a differential equation is rendered into a simpler algebraic 
expression which, in turn, becomes a fundamental equation of 
a system's functional blocks. Each block indicates the 
'relationship between its input driving signal and output 
response; ie . , the transfer function of the input and 
output. The transfer function is defined as the ratio of 
the Laplace transform of the output to the transform of the 
input with all initial conditions assumed to be zero. A 
block's transfer function does not include any information 
about the internal structure of that part of a system, 
merely the transformation of a signal between two points in 
a system. Furthermore, it should be borne in mind that all 
physical systems have certain transfer characteristics that 
are actually nonlinear to some extent. A dc motor system 
can be fairly accurately described using the transfer func- 
tion approach to modeling and, when nonlinear transfer char- 
acteristics exist, CSMP modeling can usually be obtained 
from preprogrammed nonlinear function routines. 

C. DEVELOPMENT OF MOTOR SYSTEM EQUATIONS 

In order to construct the block diagram of a system, 
one approach is to identify the various mathematical rela- 
tionships between the components of the system and to write 
balance equations that will identify each individual block 
comprising the system. The system equations for a brushless 
dc motor are presented in chapter three where all the 
details of the logic control and switching of transistors 
for commutation of the motor are developed. What follows is 
a simplified input-output characteristics model of the type 



that has been a standard for brush-type Jc motors used in 
control system studies [Ref. 3]. The development of this 
input-output model is the first learning step in the 
modeling process for the brushless dc motor and has the 
additional usefulness of being a tool for a concurrent 
mechanical engineering study cf load torque requirements for 
the motor [ Ref . 1 ]. 

To write the electrical balance equation describing a 
basic dc motor, the classical loop method based on 
Kirchoff's voltage law is applied to the equivalent motor 
circuit (Figure 2.1). 



R 




Figure 2. 1 Equivalent DC Motor Circuit 



The basic dc motor has the coil windings on the armature and 
the magnetic field source on the stator. In the case of the 
brushless dc motor/ the rotor consists of permanent magnets 
that produce the magnetic field and the stator contains the 
coil windings. The Kirchoff Law for the stator circuit is 



es (t) = Ldi/dt + Ri (t) + eb(t) (egn 2.1) 

where R and L are the stator resistance and inductance 
respectively, and eb (t) is the back emf. When a conductor 
moves in a magnetic field, or there is any relative motion 
between the magnetic field source and the conductor, a 
voltage is generated across the terminals of the conductor. 
In the case of the dc motor, the voltage is proportional to 
the shaft velocity and tends to oppose the current flow. 
The relationship between the back emf and the shaft velocity 
is 



es (t) = *wm ( t) = Kb*wm (t) 



(ey-n 2.2) 



es (t) = Kb*d9m (t) /dt 



(egn 2.3) 



Substituting and rearranging equations 2.1 and 2.2 into 
their differential forms gives the following two equations. 



di/dt= (1/L) *es (t) - (E/L) *i (t) - ( 1/L) *Kb*vm (t) (eqn 2.4) 



dem(t)/dt = wm(t) (eqn 2.5) 

The basic balance equation that describes the motor as 
a mechanical system derives from application of Newton’s 
laws of motion to the system components. The dynamic equa- 
tion for a motor coupled to its load is 

tm (t) = J*dwm(t)/dt + 3*wra (t) + tl(t) (eqn 2.6) 
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where tl (t) is the load torque. Rearranging the equation it 
hecomes 



dwm (t) /dt= (1/J) *tm (t) — (1/J) *tl (t) - (B/J) *wm (t) (egn 2.7) 

E is the total viscous friction coefficient of the motor and 
is comprised of the sum of the viscous friction due tc the 
motor. Bin, and the load, Bl, as seen by the motor shaft. 
Likewise, J is the total system inertia comprised of the sum 
of the motor inertia, Jm, and of the load inertia, Jl, 
reflected through the coupling device to the motor shaft. 
Eoth Bl and Jl are functions of the motor's speed reducing 
mechanism where N is the speed reduction ratio (N>1). The 
following equations summarize the foregoing. 

3 = Bm + Bl (eqn 2.8) 



J = J m + Jl 



{eqn 2.9) 



Bl = Blp/N 



(eqn 2.10) 



Jl = Jlp/N (eqn 2.11) 

Jlp and Blp are the inertia and viscous friction seen at the 
load side of the system. Since a dc motor is essentially a 
torque transducer that converts electrical energy into 
mechanical energy, the following equation is necessary to 

show the relationship amoung the developed torque, tm(t), 
the air gap flux, (J), and the stator current, i (t) . 

tm (t) = (Km*<J)) *i (t) (eqn 2.12) 

Because the magnetic field in the motor is assumed uniform 
and constant, this can be simplified to 
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tm ( t) = K t*i (t) 



{egn 2. 13) 



where Kt is the torgue constant of the motor. 

Eguations 2.1 through 2.13 represent the cause and 
effect eguations of the system. The application of es (t) 
produces a current flow which causes the torgue and the back 
emf tc be generated. The torgue produced then causes zhe 
angular displacement 6m (t) . Given the differential egua- 
tions of the system, the Laplace transform can be applied 
and the block diagram of the system generated. As an 
example, the block fcr eguation 2.1 is generated as follows 

es (t) - eb (t) = L* di/dt + R*i(t) (egn 2.14) 

/{es(t) - eb (t) } =X{I'* ( ii/dt + R*i(t)} (egn 2.15) 

En(s) = Es (s) - Eb (s) = (Ls + R) *1 (s) (egn 2.16) 

I(s)/En(s) = Is/ (Es - Eb) = 1/ (Ls + R) (egn 2.17) 

Performing the same operation to the remaining differential 
eguations, the block diagram in Figure 2.2 results. 

It is important to note that while a dc motor is basi- 
cally an open-loop system, the back emf of this type of dc 
motor acts as a natural feedback loop that tends to improve 
the stability of the motor. 

Cnee the basic block diagram is developed, writing a 
CSMP program is simple and fairly straightforward. Care 
must be taken, however, to put the block transfer functions 
into the proper format for CSMP's functional blocks. A copy 
of the CSMP program for a basic dc motor is provided in 
Appendix A. 
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Figure 2.2 Block Diagram of DC Motor 
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III. CURVE PRE DICTION 



A. OVERVIEW 

Cne of the purposes of this model is to produce a set 
of performance curves that accurately simulate a given 
motor's behavior when operated under varying load condi- 
tions. As such the model can be a useful tool for studying 
the changing performance due to the configuration changes of 
a motor under development or for observing the behavior of a 
motor under varying conditions of application or environ- 
ment. Thus, having developed an input-output computer model 
for the DC motor, the next step was to assign values to the 
model parameters and run the program. Program inputs, 
constants and parameters were assigned values that are 
typical of brushless DC motors commercially available. 
Table I provides the parameter values for a given motor that 
will be used as a generic motor for the following analysis 
and discussion. 



TABLE I 

Typical Commercial Motor Parameters 



Stator Resistance, R 
Stator Inductance. L 
Torque constant, Kt 
Back EMF constant, Kb 
Rotor Inertia, Jm 
Viscous Friction 

Coefficient, Em 



2 . 74 
0.0016 
15.9 
0 . 112 
0 . 001 



ohms 

henries 

oz-in/amp 

volt/rad/s 

oz-in/s 2 



unknown oz-in/rad/s 
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B. SPEED VS TORQUE CURVE 



Specifically, three curves were produced: motor speed 

vs. generated torque, current vs. torque and output power 
vs. torque. Since permanent magnet DC motors have linear 
characteristics, a straight-line curve resulted under all 
load conditions when the parameters of Table I were inserted 
into the model and a voltage of 30 VDC was applied. In 
order to fully understand the individual contributions of 
the various motor parameters on the speed-torque curve, it 
was necessary to determine the mathematical r elaticnships 
between the parameters. Referring to equations 2.1 and 2.5, 
the electrical and motor equations when considered under 
steady state conditions become 

es (t) = R*i(t) + Kb*wm(t) (eqn 3.1) 

tm (t) = Kt*i (t) (eqn 3.2) 

Solving for i(t) and combining the equations results in 

es (t) = (R/Kt) *tm (t) + Kb*wm(t) (eqn 3.3) 

When there is no torque applied to the motor, the current is 
small enough to neglect and the no-load velocity becomes 

wnl = es/Kb (eqn 3.4) 

Given that the applied voltage is a constant, it is apparent 
that knowledge of the value of the back emf constant, Kb, 
will provide the various no-load speeds for any value of 
applied voltage. On the other hand, in the absence of any 
known motor specifications, the no-load speed of a given 
motor can be empirically measured and then the back emf 
constant, Kb, can be mathematically derived. 
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looking next at tte slope of the speed-torgue curve, it 
can be seen that whereas one endpoint of the curve termi- 
nates at no-load speed, the other ends at stall which is 
that value of torgue which is large enough to literally stop 
the motor. From eguations 3.1 and 3.2, with wm set to zero 
the generated torgue, tm(t), at stall is ts given by 

ts = (Kt/R)*es (egn 3.5) 

Combining eguations 3.4 and 3.5, the eguation of a straight 
line for the speed-tcrgue curve results in 

wm (t) = wnl - (Rm*tm(t)) {egn 3.6) 

wher e 



Rm = R/(Kb*Kt) (egn 3.7) 

and Rm is the slope factor for the curve, also called the 
speed regulation constant. Obviously, changing any or all 
of the parameters will alter the slope factor, Rm. However, 
the back enf constant, Kb, and the torgue constant, Kt, are 
strictly related by the common factor, air gap flux, 4>. In 
fact, when both constants are put in the MKS system of 
units, Kb in volt/rad/s and Kt in Nm/amp, then they are 
egual as eguation 3.8 shows. 

Kb ( vol t/rad/s) = Kt (Nm/amp) (egn 3.S) 

Thus, knowledge of either constant provides knowledge of the 
other. Referring to eguation 3.7, it can be seen that the 
final unknowns are the resistance, R, and the slope factor, 
Rm. If the resistance, R, is known from manufacturer speci- 
fication sheets, for example, then the slope, Rm, can be 
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easily solved. On the other hand, if given an empirically 
measured speed-torgue curve, it is now obvious that the 
curve not only provides Kb which in turn provides Kt, but it 
also can provide a derived value for motor resistance, R. 
Cnee acquired, the three basic parameters, Kb, Kt and R, 
were plugged into the CSMP model and the speed-torgue curve 
f or • a 30 VDC input was generated (see Figure 3.1). 

This same procedure can be used to produce a family of 
speed-torgue curves if it is desired to study a motor's 
performance under different input voltages. Figure 3.2 
shows the family of curves produced when the values given in 
Table I were entered into the model. 

C. CURRENT VS TORQUE CURVE 

With the speed-tergue curve understood, the next task 
was to repeat the same general procedure for the motor 
current vs. torque. As a starting point, the equation for 
the steady-state current is given by 

iss (t) = (es(t) - Kb*wm(t))/E (egn 3.9) 

Since Kb and R were fixed by the speed-torgue curve, gener- 
ating a current-tor gue curve was straightforward procedure 
for any applied voltage (see Figure 3.3). 

D. OUTPUT POWER VS TCRQOE CURVE 

The final performance curve to be studied was the 
output power vs. torgue curve. The basic equation for power 
is given by 

P (t) = tm(t)*wm(t) (egn 3.10) 

P(t) = Kt*i(t) *wm (t) (egn 3.11) 
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Figure 3.1 Motor Speed vs Torque Curve 



where power is measured in watts. Since this equation 
involves no new constants, producing the power curve was 
only a matter of adding a line to the existing model which 
also included a conversion factor for changing the torque 
units from oz-in/amp to Nm/amp in order for the power to be 
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Figure 3.2 Family of Speed-Torque Curves 

dimensionally correct in watts. Again running the program 
with the parameters given in Table I, the model produced the 
power curve given in Figure 3.4 
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Figure 3.3 Motor Current vs Torque Curve 

As expected both the speed-torque and curren t- torque 
curves were linear. The power-torque curve is obviously not 
linear, especially at high leads where the power begins to 
roll off fairly sharply. Fro® the modelling point of view 
this is due solely tc the fact that the motor slews down 
faster than the load torque increases (see equation 3.9) . 
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OUTPUT POWER VS TORQUE CURVE 




Figure 3.4 Output Power vs Torque Curve 

Another factor that further accentuates this power rolloff 
in an actual motor is the effect due to armature reaction. 
Whenever a current flows through the motor armature in a 
permanent magnet dc motor/ the armature becomes an electro- 
magnet which produces a flux that tends to oppose the flux 
produced by the permanent magnets. This has the effect of 
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partially demagnetizing the permanent magnets. This 
demagnetization is a reversible effect, meaning that as 
current returns to zero the permanent magnets return to 
their full strength. However, as the current increases and 
armature reaction sets in, it has the effect of decreasing 
the torque constant, Kt (as well as Kb). As equation 3.9 
indicates, this also contributes to the output power rolicff 
at high loads. Though current permanent magnet DC actors 
made of rare earth magnets have high coercivity that resists 
armature reaction, at high levels of rated current a small 
amount of armature reaction still occurs. For the purposes 
of this thesis, the model is considered to accurately simu- 
late motor behavior fcr loads up to the peak power load tor 
the given applied vcltage. From that load upward, it is 
assumed that armature reaction is likely to occur and result 
in nonlinear behavior which is not included in this model. 

E. MCTOB REVERSAL 

The final step in the modeling procedure was to reverse 
the motor’s direction and ensure that mirror images of the 
three curves resulted. To do this, it was first necessary 
to replace the positive input voltage with a negative count- 
erpart. Next, in order to maintain model consistency where 
the lead terque, tl(t), is treated as a torque that opposes 
the motion of the motor, it was necessary to modify the 
model so that the load torque was a positive value. This 
consisted of creating a C3HP procedure block that identified 
the input voltage as either positive or negative and then 
treated the load torgue accordingly. Once program tugs were 
removed, the model accurately simulated motor reversal with 
proper speed, current and power values for the appropriate 
load conditions. 
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IV. ELECTRONIC COMMUTATION 



SWITCHING ACTION 



The next major step in this thesis was to convert the 
basic model of a standard brush-type do motor into a three- 
phase, fcur-pole brushless dc motor system. Figure 4.1 
shows this system with the three-phase stator windings 
configured in a standard ’star' connection with each winding 
oriented 120 electrical degrees from the others. 




Figure 4. 1 Controller Configuration 

The six transistors are connected to the ends of each stator 
leg and through logic-controlled switching action provide 
three-phase full-wave motor control. Figure 4.2 indicates 
the switching logic sequence for counterclockwise rotation. 

To produce clockwise rotation, the sequence is reversed. 
In either case, a pair of transistors will be switched at 
the start of each 30 degree (mechanical) interval which 
causes current to simultaneously build up in one leg, flow 
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Figure 4.2 Switching Logic of a 3-Phase Brushless DC Motor 

steadily in a second, and decay to zero in a third. Given 
proper sequencing, the net developed torque ideally readies 
a steady state value that causes motor rotation at a 
constant speed in the desired direction. 

For windings configured in a star arrangement, it can 
be seen that conduction is continuous in one leg while 
commutation occurs in the other two legs. For example, as 
Q 1 and Q5 are energized between 30 and 60 degrees (mechan- 
ical) , current flows down leg A and into leg B. At the same 
time, the current through leg C due to the energy storage 
behavior of an inductor decays to zero through diode, D3 
(see Figure 4.1). In the next seguence, 60 to 90 degrees, 
Q6 is energized and C5 switched off. Current has reached a 
steady state flow through leg A, but now it flows down leg C 
and the current in leg B decays to zero through D2. In the 
next sequence, 60 to 120 degrees, Q2 is switched on and Q1 
switched off. Current has reached steady state flow through 
leg C but now current builds up and flows from leg B while 
the current in leg A decays through leg C. This same action 
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repeats in subsequent commutation states. Looking at point 
E in Figure 4.1, it can be seen that, since current is 
always continuous in one leg whether flowing into cr from 
node D, Kirchoff's Current Law requires that the net current 
flowing in the other two legs must equal the flow in the 
continuous leg. This permits the following important 
simplification which was needed later in the modeling 
process. Namely, the three-legged stator can be approxi- 
mated by just two windings serially configured. One leg 
will always accurately reflect steady state flow conditions 
while the other will be treated as if it were in steady 
state because the build up in one leg is balanced by decay 
in the other. It is recognized that this model is a first 
approximation and that future modeles will require addi- 
tional refinements in order to adequately represent the 
effect of switching transients on the transistor and diode 
elements. 

E. HAIL EFFECT SENSCB FEEDBACK 

As was mentioned in the introduction, the switching 
action of the commutation circuitry is based upcn rotor 
position feedback from Hall effect sensors. In crder to 
represent electronic switching and postional feedback, as 
well as the three phase winding configuration, it was neces- 
sary to modify the block diagram of Figure 2. 1 to that in 
Figure 4.3 

Comparing Figures 2.1 and 4.3, it is seen that the 
chief difference is in the switching logic block and the 
addition of two more transfer function blocks to account for 
the additional windings. To better understand the logic 
relationship between the sensors and the current switches, 
the truth table for both counterclockwise and clockwise 
rotation is presented in Table II. 
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I 

Figure 4.3 Equivalent Circuit With Electronic Commutation. 



A brief review of the action af a Hall effect sensor is 
in order at this point. Consider a conducting material with 
a constant current flow. When there is no external magnetic 
field present, the current flows undisturbed in a straight 
line. But when a magnetic field is applied, a deflecting 
force, the Lorentz force given by II x B (where I is the 
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TABLE II 

Sensor and Switching Logic 



BPS A 

1 

1 

1 

0 

0 

0 



RPS B 

1 

0 

0 

0 

1 

1 



BPS C 

0 

0 

1 

1 

1 

0 



PH A 

high 

high 

open 

low 

low 

open 



PH B 

low 
open 
hi gh 
high 
open 
lo w 



PH C 

open 

low 

1 cw 

open 

high 

high 



COUNTERCLOCKWISE 



0 

1 

1 

1 

0 

0 



0 

0 

0 

1 

1 

1 



low 


high 


open 


open 


high 


1 cw 


high 


open 


low 


high 


low 


open 


open 


low 


high 


low 


open 


h i c h 



CLOCKWISE 



current, 1 is the length ox the conductor and 3 is the 
magnetic field) causes the current to bend and electrons to 
pile up on one side and positive charges, or holes, to do 
the same on the other side. Thus a transverse Hall poten- 
tial difference, Vh, develops - across the conductor (see 
Figure 4.4). 

Figure 4.5 shows that in the presence of a sinusoidal 
magnetic field produced by a rotating magnet, such as one of 
the pole-pairs of a multiple-pole motor, the induced Hall 
voltage takes on an ideal square wave shape. 

The logic circuitry of the sensor is configured in such 
a way as to output a logic level 1 for positive Hall volt- 
ages and a 0 for negative voltages. Briefly reviewing the 
action of the motor: as the shaft turns, the rotor position 
sensors send this information to the switching logic block 
which then decides which phases are to be made high, lew or 
open in accordance with Table II. Equivalently, the 
voltage, En, is then applied to the two appropriate windings 
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Figure 4.4 The Hall Effect 

and zero voltage is applied to the third winding in agree- 
ment with the assumption that current is steady in one leg 
while the commutation in the other two legs is eguivalent to 
a continuous current in one leg. As before, each voltage 
then produces a current. Because this is assumed to be a 



33 





F.PS 



1 




Figure 4.5 Sensor Logic Based on a Hall Effect Sensor 

linear system, the principle of superposition is applicable 
and, thus, the torques developed by current flow through 
separate windings car be considered independently and then 
summed to give a total motor torque. Beyond this point, the 
brushless motor components are identical to its brush-type 
counterpart and behave similarly. 
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C. KCDEI REVISION 

CSilP modeling of this version of the brushless do motor 
called for two modifications to the brush-type model. 
Essentially, a CSMI procedure function was constructed to 
act as a subroutine that embodies the logic of the switching 
action. The logic recognizes that if the motor shaft is 
within seme degree range, for example 30 to 60 degrees, then 
it sets the appropriate rotor position sensors and energizes 
the switches for the appropriate two windings and deener- 
gizes the switch for the third. In effect, only the forcing 
voltages are passed to the main program, whereas switch and 
sensor positions are available for output (see Appendix A). 
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V. VARIABLE FLOX 



A. AIR GAP FLOX 

The steady state behavior of both a brush-type and a 
brushless dc motor has been studied under the fundamental 
assumption that the back emf generated is a product of tae 
shaft velocity times a constant, Kb (see Equation 2.2). 
Likewise, the motor torque is equated to the product of the 
current times a constant, Kt (see Equation 2.12). But as 

Equation 2.11 reveals in the case of the torque constant, 

Kt, the torque is in fact the product of some other constant 
cf proportionality, Ktl, times the air gap flux, <j). Ihe 
same is true for the back emf constant, Kb; that is, Kb is 
the product of some proportionality constant, Kbl, times the 
air gap flux, <J). In both cases, the flux, $, has been 

assumed to be constant which is in agreement with the 

general practice when . modeling permanent magnet dc motors. 
In point of fact, the flux is not constant but rather is 
distributed sinusoidally, generally being the sum of a sinu- 
soid and one or more of its harmonics. How the motor is 
physically constructed determines this flux distribution. 
Two of the chief factors are the structure of the magnetic 
pole faces and, most importantly, how the windings are 
distributed in the armature. The latter includes such 
factors as the number of slots, the number of coils psr slot 
and the spacing of the slots [Ref. 5]. The careful designer 
distributes the windings so as to minimize the effects of 
harmonics. Regardless of their constituents, these 
composite sinusoids have an average value, of course, and so 
in practice it is this average value that becomes the 
constant flux term, <t • 
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B. FLUX AS AN AVERAGE VALUE 



A logical next step in the design of a more realistic 
model was the determination of an average value fur the 
varible flux, ({)avg, and the corresponding proportionality 
constant, KK1 such that the ratio of <()avg and KX 1 equals cae 
and the relationships of equations 2.2 and 2.5 are 
preserved. The equations for back emf, eb (t) , and developed 
torque, tra (t) , now become 



eb (t) = Kbp * wm (t) 


(e qn 5.1) 


tm (t) = Ktp * i (t) 


(eqn 5.2) 



where the back emf constant, Kbp, and the torque constant, 
Ktp, are now expressed as 



Kbp = Kb * (<{>avg/KKl) 


(eqn 5.3) 


Ktp = Kt * ((J)a vg/KK 1) 


(eqn 5.4) 



Assuming as a first approximation that the air gap flux 
varies as a simple sinusoid of unit magnitude, the first 
step toward finding the total system flux, $avg, was to 
determine the relationship of the individual flux patterns 
in each of the three phase windings. Analysis of a four- 
pole magnet rotating under three windings separated by 120 
degrees (electrical) revealed the following relationships 



<j)a = sin (26 + tt/6) 


(eqn 5.5) 


(J)b = sin (20 + 91V6) 


(eqn 5.6) 


<])c = sin (2© + 5T1X6) 


(eqn 5.7) 
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where <t>a, <J)o and 4)c are the flux contributions of the 
respective phases. returning to an earlier assumption that 
the circuit can be approximated as two windings in series, 
ignoring the third, then the average flux can also be 
approximated from the sum of the flux developed in two wind- 
ings at a time. This approximation is further supported by 
the fact that the flux in the decay leg has, in fact, an 
average value of zero during every commutation interval. 
From the switching sequence of Figure 4.1, the two 
conducting legs were identified for each interval and, from 
Equations 5.5 through 5.7, their respective flux values 
calculated and summed. For every interval, the resulting 
flux waveshape was identical and had an average value of 
(Jiavg = 1.631 (see Figure 5.1). 

Since KK1 merely normalizes (|)avg, its value was easily 
determined to be KK1 = 0.613. The appropriate substitutions 
were then made in the CSMP model (see Appendix A for 
Revision Two). Since this amounted to nothing more than a 
rearrangement of constants, the model's output behavior 
remained unchanged. At this stage, the model simulates the 
behavior of a brushless dc motor being electronically commu- 
tated and exhibiting constant, steady state developed torque 
and motor speed due to a fixed, average value air gap flux. 

C. SINUSOIDAL FLUX 

The final stage of this thesis was to examine the 
steady state behavior of the brushless dc motor under vari- 
able air gap flux conditions. As a first approximation, the 
flux was modeled as a simple sinusoid and the flux relation- 
ships are those given in equations 5.5 through 5.7 As 
before, only the flux relating to the two conducting legs 
was considered. This composite variable flux, <J)v 1 , is shewn 
in Figure 5.1 Equations 5.3 and 5.4 had to be modified to 
account for the different value of flux. 
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FLUX VARIATION 




Figure 5.1 Composite Flux Variation for Two Windings 

Kbpp = Kb * (<J)v1/KK2) (egn 5.8) 

Ktpp = Kt * (Cf>v1/KK2) (egn 5.9) 

Recall that the objective was to keep the values of Kbpp and 
Ktpp as near as possible to those of Kb and Kt, respec- 
tively. In order to roughly normalize (f) v 1 , the constant 
value was determined as KK2 = 1.655. Since ({)v 1 was now a 
variable quantity, it should be recognized that Kbpp and 
Ktpp were no longer strictly constants and were better 
expressed as the following functions. 

Kbpp = K b (4>v 1 ) (egn 5.10) 
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Kt pp = Kt (<f)v 1) 



( egn 5. 11) 



The equations for back emf , eb (t) , and developed torque, 
tm (t) , now became 



eb (t, <J) v 1 ) = 


Kbpp 


* wra (t) 






(eqn 5. 12) 


tm (t, <f> v 1 ) = 


Ktpp 


* i(t) 






(egn 5. 13) 


Eecause 


the 


total flux 


was no 


longer 


a constant. 


average value. 


it 


was expected 


that 


variation 


in developed 


torque (and. 


consequently, in 


motor 


speed) 


would result. 



Furthermore, as the system block diagram of Figure 2. 1 indi- 
cates, torque is converted by the motor transfer function 
into a rotational speed. This is a process that is equiva- 
lent to inertia filtering of the torque ripple that results 
from variable flux. For the purposes of this thesis, ripple 
is defined as the ratio of the amount of variation from the 
maximum to the maximum itself. Having mace the appropriate 
substitutions in the model (see Appendix A for Revision 
Three) , the program was run under .no-load conditions. It 
was observed that there was, in fact, ripple in both the 
developed torque and motor speed and that some inertia 
filtering did occur. The next logical step was to lead the 
motor and study the effect on torque and speed ripple. 
Table III indicates that under no-load conditions where 
speed was highest, the percentage of ripple in the output 
speed was at a minimum but that torque ripple was at a 
maximum. 

It can be seen that as the load increased the amount of 
ripple in the developed torque decreased from about 8 4T to 
8%. A closer examination of the torque behavior under 
varying loads revealed that, even as the motor torque neces- 
sarily grew with increasing loads, in every case there was 
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TABLE III 

Torque and Speed Ripple Due to Sinusoidal Flux 



MOTOR 


MOTOR 


% RIPP 


LE 


LOAD (oz-in) 


SPEED (RPM) 


TORQUE 


SPEED 


0.0 


2557 


33.9 


1. 1 


32.0 


2087 


2 1. 1 


1.6 


€«. 0 


1617 


12.2 


2.6 


96.0 


1147 


7.8 


4.6 



approximately a constant 8-9 oz-in variation. This varia- 
tion can be attributed primarily to the flux variation in 
Ktpp and to a lesser degree to the variation in i(t) caused 
by the variation in back emf. Thus the magnitude of torque 
variation remained about the same for all loads though the 
percentage ripple appears to indicate otherwise. In the 
case of motor speed, the amount of ripple appeared to 
increase slightly with increasing load and thus decreasing 
speed. In this instance, the magnitude of variation did 
increase a small amount with decreasing speed - from 28 rpir. 
at no-load to 52 rpm at 96 oz-ins - and suggests that speed 
ripple, at least, is somewhat dependent on load. This is 
due in part to the fact that as the motor slows, there is 
less rotor momentum and thus less inertia filtering. In 
addition, at slower speeds there are less frictional effects 
which at higher speeds also tend to filter out ripple. In 
general, it makes intuitive sense that at higher speeds 
frequency variations are harder to discern than at slower 
speeds. 

D. HARHCNIC FLUX 

Having observed the effects of simple sinusoidal flux 
variation, the next step was to replace the single sinusoid 
with a flux composed of the sum of a sinusoid and its 
harmonics which is mere representative of the air gap flux 



patterns in an actual motor. The back emf of a typical 
commercial brushless dc motor is given in Figure 5.2 where 
the voltage between two terminals was measured with the 
motor rotating at 1200 rpm. 

A harmonics analysis of this waveshape shows that the 
principle harmonic was the fifth and that the waveshape was 



closely described by the expression 




3.0 sin (0) + 0.59 sin(56) 


(eqn 5. 14) 


Ehen substituted into equations 5.5 thru 5. 


,7, and the 


different angular speed taken into account, 
relationships for flux result 


the fc Hewing 


4) a = 3. 0*sin (20+ tr/6) + 0. 59*sin ( 1 00 + 5 ir/6 ) 


(eqn 5. 15) 


4>b = 3.0*sin (28 + 9 7T/6) + 0 . 59 *sin { 1 09 + 9 fl/6) 


(eqn 5. 16) 


<|>c = 3.0*sin {26+5 TT/6) + 0 . 59*sin ( 1 06+ 77/6) 


(eqn 5.17) 



The same procedures used in the simple sinusoid model 
were then applied to the harmonics case where the composite 
flux changes to (J>v2, a flux factor that varied to a slightly 
greater degree than did <{>v1. Equations 5.3 anu 5.9 had to 



be modified as follows 




Kbppp = Kb * (<j)v 2/KK3) 


(eqn 5. 13) 


Ktppp = Kt * ((pV2/KK3) 


(eqn 5. 19) 


where KK3 approximately normalizes the flux 


and had the 



value KK3 = 5.38. As before, Kbppp and Ktppp varied with 
the flux and so the new equations for back emf and developed 
torque became 
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MAGNITUDE 



BRCK EM r WAVEFORM 



o 




Figure 5.2 Back EMF Waveform Due to Harmonic Flux 
eb(t,<J)v2) = Kb ($v2) * wm(t) (egn 5.20) 



eb(t,(fv2) = Kbppp * win (t) 



(egn 5.21) 



tm (t,(J)v2) = Kt ((J) v 2 ) * i (t) 



(eqn 5. 22) 



tm(t,(f)v2) = Ktppp * i (t) (egn 5.23) 

The appropriate substitutions were once more made in the 
model (see Appendix A for Revision Four) and the model was 
run under various loads. As anticipated, the variable flux, 
<fv2, produced ripple in the developed torque and motor 
speed. Table IV shews the resultant values for various 
loads . 



TABLE IV 

Torque and Speed Ripple Due to Harmonic Flux 



MOTOR 


MOTOR 


% 


RIPPLE 


LOAD (oz-in) 


SPEED (rpm) 


TORQUE 


SPEED 


0 . 0 


2557 


94.6 


2. 4 


32.0 


2087 


33.5 


3. 7 


64. C 


1617 


22.5 


5.6 


96.0 


1147 


14.7 


9.8 



Comparing Tables III and IV, it can be seen that the values 
for torque and speed ripple are of similar magnitudes. The 
somewhat higher values for torque are due to the fact that 
in the <fv2 case there was a 16-21 Ox-in variation as opposed 
to an average 8-9 cz-in variation in the <{)v1 case. The 
speed ripple behavior was much like the (|) v 1 case with 
proportionate increases resulting from the increased torque 
ripple. Thus, for both torque and speed ripple due to <J>v2, 
the model performed as expected; that is, slightly higher 
values of ripple occurred which were clearly attributable to 
the increased variability of the flux pattern of composite 
sinusoidal harmonics. 

An important observation made during runs of the <{)v1 
model came as a consequence of the motor speed variability. 
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Since no-load 



particularly when the motor was unloaded, 
speed was a function of Kb as given in equation 3.4, the 
variable Ktpp had the effect of increasing the upper-end 
speed beyond the 2557 rpm baseline given in Figure 3.1 This 
increased magnitude ccupled with the swing of speed values 
created an undesirable effect in the unloaded condition. In 
effect, at certain points in time the additional speed 
produced bach emf values that exceeded the input voltage 
which in turn resulted in negative current and developed 
torque values. Though these were transient negative values 
and had small effect on the steady state operation of the 
motor, still their presence was obviously unrealistic. The 
solution to returning the motor to all positive values of 
current and torque was to increase the friction filtering 
effect by adjustment of the viscous friction coefficient, 
Bm. It will be recalled from Chapter Three that a value for 
the Bm of the commercial motor being modelled was not avail- 
able and that for the basic model a value for Bm was derived 
from the curve fitting process. Thus, it seemed reasonable 
and permissable to repeat the procedure here. 

Since it makes sense that increasing the friction within 
a system will slew the system down, in this case increasing 
the value of Bm did just that. By increasing Em from 
0.00015 to 0.045 oz-in/rad/s, the upper-end of the motor 
speed was reduced to a value that eliminated the undesirable 
negative values. This caused the no-load speed to be 
lowered slightly below the desired value of 2557 rpms. 
Recalling that no-lcad speed is fixed by the tack emf 
constant, Kb, this parameter was adjusted from 0.1120 to 
0.1089 vclt/rad/s to bring the back up to speed. Since this 
adjustment was less than 3% and the manufacturer specifica- 
tions allowed a 10 5? +/- measurement error margin, this 
adjustment to Kb was permissable. Of course, a propor- 
tionate adjustment to the torque constant, Kt, was also made 
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since Kt and Kb must maintain a constant relationship as 
explained in Chapter Three. The effect on motor speed 
ripple was neglible. Reviewing the $v2 model at this point 
revealed that similar behavior in the no-load state occurred 
but on a larger scale due to the greater variability of §v2. 
A similar set of adjustments was necessary to eliminate the 
negative current and torgue values in the unloaded condi- 
tion. It required an adjustment of Bm from 0.00015 to 0.045 
oz-in/rad/s in order to ensure that all values were posi- 
tive. In addition, the back emf constant, Kb, was adjusted 
to 0.1192 volt/rad/s and a proportionate adjustment to Kt 
was made. Again, all adjustments were within tolerances. 

Of course, it was earlier established in Equation 3.7 
that changes to Kb and Kt would change the slope of the 
speed-torgue curve. Thus, in both cases the slope was 
returned to its original value by adjusting the motor resis- 
tance, R. In the case of sinusoidal variable flux, R was 
reduced from 2.74 to 2.70 ohms. For harmonic variable flux, 
R was increased from 2.74 to 3.0 ohms. Again, both adjust- 
ments were within the 10% +/- allowable deviation. It was 
also observed that since both Kb and Kt no longer remained 
constant but rather varied within small ranges, the slope 
was not constant and varied accordingly in both the ijivl and 
$v2 models with the greatest change occurring when Kb and Kt 
simultaneously reached the minimum values of their respec- 
tive ranges. Bearing in mind that these deviations in slope 
are fairly transient, their average values are most impor- 
tant and are the approximate values from which the final 
adjustments to R were based. 
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VI. SUMMARY 

A. REMARKS AND C0NC1DSI0NS 

CSMP language is a convenient tool for modeling a 
brushless DC motor. Positional sensor feedback and 
switching logic being functions of time, the action cf elec- 
tronic commutation can be nicely simulated in CSMF. The 
model can be fairly easily modified in order to study ether 
motor configurations. For example, the three-phase star 
configuration used in this thesis could easily be changed to 
a grounded neutral, point D in figure 4.1, and the same 
types of analysis performed. In future studies where the 
complexity of the mcdel is increased in order to include 
more design detail, the advantage of CSMP will become even 
more apparent. 

In the course of simulating the effects of different 
air gap flux variations, it appeared necessary to make minor 
parameter adjustments in order to offset the feedback 
effects cf torque and speed ripple. Within the scope of 
this thesis, adjusting the steady state performance was 
reasonable and justifiable. In a larger context, however, 
various control systems, particularly torque generation 
controllers, are specifically designed to reduce torque 
ripple and its effects. So, in at least one sense, these 
adjustments were somewhat artificial. To the extent that 
parameter adjustment is necessary and desirable, it must be 
understood that this is not a simple procedure since adjust- 
ment of one parameter almost always requires adjustment of 
one or more others. These adjustments in turn necessitate 
readjustment of the original parameter. The process is thus 
an iterative one. This model permits this kind of parameter 
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adjustment, but it is a long and tedious process that would 
he better performed ty some sort of optimization algorithm. 

E. RZCCHHENDATIONS FCR FUTURE STUDY 

The simulation cf ripple effects is particulary rele- 
vant to the control engineer who must design the control 
circuitry and power conditioner to meet specific performance 
requirements. This points to several areas for future 
study. One area is the control of torque generation and the 
several principle configurations used in brushless DC 
motors. Two commonly used methods are the sinusoidal torque 
generation principle and the trapezoidal torque scheme. 
[Ref. 5]. In addition to torque generation control, a 
second area for investigation is power control of brushless 
DC motors. In this thesis, output power was unregulated and 
was studied simply as a function of a constant input 
voltage. In practice, however, power is produced and 
controlled by varying the supply voltage. Various schemes 
exist for power control; amoung these are linear transistor 
control and pulse-width or pulse-frequency control. 
[Ref. 5]. 

A very important area that this thesis did not address 
was that of the switching transients that occur as a conse- 
quence of commutation. During each switching interval, each 
winding is an inductor that stores energy which causes a 
significant amount of current to flow at the instant of 
switching. This can have serious effects, especially if 
breakdown conditions exist amoung the commutation transis- 
tors and other control circuit elements. In this same area, 
switching transients also occur within the control circuitry 
itself. Thus, a thorough investigation and understanding of 
the behavior of power transistors, diodes and other control 
elements used with electronic commutation needs to be dene. 
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APPENDIX A 

LISUNG OF MODEL PROGRAMS 
A. BASIC PROTOTYPE MODEL 

The program in this section is a basic program that 
simulates the action of a standard brush-type DC motor. It 
is the prototype from which the brushless DC motor model and 
its several revisions are derived. 
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TERMINAL 

TITLE BASIC CC MOTOR SYSTEM 

TIMER FINTIM = .04, CUTCEL = 0.0004, PR CEL = 0.C004 
PRINT VEMF , IM, TM , NM , NMRPM, PwR 
PACE MERGE 



E. REVISION ONE 



This is the first revision of the basic prototype. This 
model simulates the switching logic of a brushless EC motor 
which is based on feedback from Hall effect sensors. In 
addition, the windings are treated independently and their 
contributions to developed torque are superposed. Since the 
superposition results in twice as much current flow as a 
single lumped coil, the total current is halved to retain 
the same overall motor behavior as in the prototype. 
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C. REVISION TWO 

In the following program, the motor’s air gap flux is 
treated as the average value of the sum of the flux acting 
in two windings as explained in Chapter Five. For conven- 
ience of study, both windings are again lumped into a single 
winding as in the prototype. This facilitates a closer 
focussing on the effects of variable flux on motor behavior. 
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D. REVISION THREE 

The purpose of this version of the program is trea 
air cap flux as varying in simple sinusoidal fashion, 
total system flux is treated as the algebraic sum o 
flux developed in the two active windings. 
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This version of the model treats the flux as varying 
according to the sum of a sinusoid of fundamental frequency 
and its fifth harmonic as explained in Chapter Five- The 
total flux is again approximated as the algebraic sum of the 
flux developed in two windings at a time. 
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THE FOLLOW IKG PROCEDURE hAS THE SAFE FCRM A KD EFFECT AS THE 
PREVIOUS. 17 IS USED FERE TO CETERMINE THE ALGEBRAIC SUM CF 
THE VARIABLE FLLX CCMPCNEMTS FCR USE IK COMPUTING THE MCTOR'S 
APPROXIMATE EACK EMF IEEMFTI. THIS FUNCTION COULD AS EASILY 
HAVE bEEN C/LCCLATEC IK TFE PREVIOUS PRCCECLRE EXCEPT THAT IT 
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APPENDIX B 



SAMPLE OUTPUT 



The following is a sample of the type of information 
the model can provide. In this particular case, the fourth 
revision of the basic prototype, the harmonic flux case, was 
used. For completeness, the motor was run at loads from 0.0 
up through 96.0 oz-in, which was the range of leads with 
which all analysis was performed. In addition to such motor 
variables as speed, current and power, a sampling of the 
feedback from Hall effect sensors as well as the switching 
action of the power transistors is included. 
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